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The damped harmonic oscillator is a workhorse for the study of dissipation in quantum mechanics. 
However, despite its simplicity, this system has given rise to some approximations whose validity 
and relation to more refined descriptions deserve a thorough investigation. In this work, we apply 
a method that allows us to diagonalize exactly the dissipative Hamiltonians that are frequently 
adopted in the literature. Using this method we derive the conditions of validity of the rotating-wave 
approximation (RWA) and show how this approximate description relates to more general ones. We 
also show that the existence of dissipative coherent states is intimately related to the RWA. Finally, 
through the evaluation of the dynamics of the damped oscillator, we notice an important property 
of the dissipative model that has not been properly accounted for in previous works; namely, the 
^q' necessity of new constraints to the application of the factorizable initial conditions. 
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I. INTRODUCTION 



, The study of dissipative systems and in particular of the Brownian motion has been pursued for a long time in 
the context of classical Q and quantum mechanics |Q. Although there has been a number of publications in this 
area there are some subtle points that have never been properly investigated in the literature. Among these we 
O^l ' could mention three major ones; a careful investigation of the relation between different models ||, the existence 
of dissipative coherent states or the condition for the employment of factorizable initial conditions. These are 
£Q ■ exactly the issues we shall address in this paper. 

Usually the dissipation in the system is described as a consequence of its coupling to a reservoir. The properties 
of this dissipative systems are generally studied through the evaluation of the time evolution of its reduced density 
operator. This evolution is often described either by a generalization of the Feynman- Vernon approach [pp|-|lC|] 
, or through master equations |i|-^,|ll|-|l5|l . In this work the properties of the system will be studied through exact 
diagonalization of different Hamiltonians of the dissipative models. 

We will consider a damped harmonic oscillator. The usual models of dissipation consist of coupling the harmonic 
oscillator to a reservoir that is conveniently chosen as a group of N noninteracting oscillators. The coupling between 
the two systems is bilinear in the creation and destruction operators of quanta of energy. Then the Hamiltonian of 
Q H . the total system is given by j| 



being 



H — Hsis + HRes + Hint, (Id) 

H S is = Tiojod^a, H Res ^ Ti'^ujjb^bj, (1.2) 

3 

H Int =h(tf + a)J2 (Wi + k ?l) , (1-3) 

3 

where we consider a harmonic oscillator with frequency lu (the system of interest) interacting with a bath of oscillators 
with frequencies ujj through the coupling constants kj's. We will take the limit of a continuous spectrum of excitations 
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in the reservoir of the Hamiltonian H . Then we will diagonalize H and determine the time evolution of the operator 
a exactly. The analysis of a (t) will determine the conditions of validity of the rotating- wave approximation (RWA) 
which consists of neglecting the terms kj&bj + k*a^b\ in (O) and writing 



H?£ A =hY^(k j & r b j + k*aV j ), (1.4) 

j 



Once this has been accomplished we will discuss the existence of dissipative coherent states. Some authors 
have stated that the coherent states are special states that remain pure during its decay in dissipative systems. We 
will show that the existence of these dissipative coherent states is directly related to the RWA; they can only exist at 
zero temperature and in systems that meet the conditions required for the RWA. 

Once we have determined the evolution of the operator a (t) of the system we can determine the evolution of any of 
its observables. However, the dynamics of these observables will depend on the specific form adopted for the coupling 
constants kj as functions of the frequencies ujj . Our method holds for an arbitrary form but in order to compare our 
results with the Caldeira-Leggett model [^J, we will reduce our results to the case where the function becomes the 
same as the one they have adopted. Then, as in Refs. [^,0, we will determine the evolution of the mean value of the 
position operator (q (t)) of the damped oscillator. The result of this calculation reveals a very special need to carefully 
treat the initial time of the motion. We propose a simple initial condition that eliminates the initial transient that 
would appear in the evolution of (q (t)) and we believe that it is enough to eliminate most of or maybe all the initial 
transients (in a certain time scale) which were noticed in this system in previous works . 



The paper is organized as follows. In Sec. II we write the Hamiltonian (LI) in the limit of a continuous spectrum 



for the reservoir exci tatio ns and we diagonalize it exactly within and without the RWA. We compare the model given 



by the Hamiltonian (IT) with the dissipative model presented in |2| in Sec. III. Here, we also determine the relation 



between the coupling function |w(oj)| 2 , introduced in Sec. II, and the spectral function J {to) introduced in |8|]. In 
Sec. IV we analyze the relevance of the different terms that appear in the calculation of the evolution of the operator 
a (t) with relation to the intensity of the dissipation in the system. In Sec. V we show under which conditions the 
evolution of the operator a (t) is reduced to that given in the RWA. In Sec. VI we show that the existence of dissipative 
coherent states is only possible within the RWA. In Sec. VII we present the calculation of the evolution of the mean 
value of the position of the damped harmonic oscillator. In Sec. VIII we discussed the physical meaning of the initial 
condition proposed in Sec. VII. Finally, we discuss the main results and conclusions in Sec. IX. 



II. DIAGONALIZATION OF THE DISSIPATIVE HAMILTONIANS 



A. Treating a Reservoir with Continuous Spectrum 

We can rewrite the Hamiltonian ([Ll]) considering a continuous spectrum of excitations in the reservoir by making 
use of the transformation between the discrete boson operators bj and the continuous ones &n 16 



dtlbn, 



(2.1) 



l/9(^) 



where g (flj) dfij is the number of modes in the reservoir with frequencies between Qj and Clj + dfij and Ji/ g ^Q .^ dft 



represents an integration in a band of width l/g(Qj) around flj. 
relation 



bnM 



5 n-n 



The operators &n then satisfy the commutation 

(2.2) 



and all other commutators v anish . 

Under the transformation (|2.lD we find 



H 



Int 



h (a f + a) dny/g(n) \k (ft) b n + k* (Q) 'b\ 



(2.3) 



where we considered that g (fij) and k (ftj) are constant inside the interval 1/g (ftj) and that Y) j fi/ g (n ) ^ ^ s nothing 

but J dfl, where this last integral covers the whole spectrum of excitations of the reservoir. Then the total Hamiltonian 
of our system is given by 



2 



H = Tiijj a) a + fi J fiS^&^dfi 



(2.4) 



+h (a T + a) 



v (fi) b n + v* (fi) b n 



dfi. 



where 



w (o) = vgljrjk (fi) 



(2.5) 



B. The Hamiltonian within the Rotating Wave Approximation 

We will now perform a canonical transformation and apply the procedure proposed by Fano [fl7f in order to 
diagonalize the Hamiltonian of our global system in the RWA that is written as 



v (fi) a)b n + v* (fi) ab 



(2.6) 



dQ. 



The diagonalization procedure presented in the sequel is basically a review of the method presented in |T^] . Our goal 
is to find an operator that satisfies the eigenoperator equation 



RWA 



Till! A,, 



(2.7) 



and therefore has its evolution trivially given by A^ (i) = A u e luJt . 

The new operator A^ can be written in terms of the operator a of the system and of the operators 6q of the reservoir 
in the form 



dfi/^n&r 



(2. 



Substituting this expression for A^ as well as (2.6) for H RW in (2.7) and calculating the commutators we have 



ujoa^a + a u J dSlv (fi) bn + J dfifi/^n&n 

+ J dQv* (fi) Pcj^a = lu [a u a + J dfiA^n^fi) 

Now, taking the commutator of this expression with and b n , we obtain 



UoOLu + J dQv* (fi) /3 Wl n = wa u , 
v (fi) + fi/3 Wl n = ufiu>,n, 



respectively. Imposing 



we have 



= S (w — Co) , 



The system of equations (2.10), (2.11) and (2.13) is identical to the one presented in 



[U! - L0 o - F (oj)] 2 + 7r\v(uj)\ 2 



(2.9) 

(2.10) 
(2.11) 

(2.12) 
(2.13) 

The solution is given by 
(2.14) 



3 



with an arbitrary phase of a u , and 



V 



1 UJ — UJ a — F (uj) 



uj — SI 



v(Sl) a u , 



where 



F(uj)=V 



\v(si)\< 



- si 



dSl, 



(2.15) 



(2.16) 



and V denotes the principal part. 

We can calculate the evolution of the operator a of the system expressing it as function of the operators A w . We 
can write a as function of A u in the following way 



(2.17) 



Taking the commutator 
operator a is given by 



a, At 



first using (2.8) and then (2.17), we obtain f u = a*. Therefore the evolution of the 
a(t) = I dioalKe- 1 ^. (2.18) 



Substituting the expression for A w in this equation and using (2.15) we obtain 

a(t) = J duj\aj\ e'^d 

1 



+ / dOw(fi) ^ / du\a u \*V 



- si 



(2.19) 



+ - 



an 



\v(m 



[n-u -F (SI)} e^ s " ^ 6 n 



C. The Hamiltonian without the Rotating Wave Approximation 



Now we will present the diagonalization of the Hamiltonian (2.4) without the RWA. The procedure that we will 
present is similar to the one adopted in [ pj[ . 

Again we want to find an operator A^ that satisfies (2.7), with H in the place of H RWA , and ( 2.12] ). Then we write 
A,,, in the form 



A u ,=a u) d + / dVL(3u,nbn + Xw<$ + / d£la u ,nbl, 



Imposing (2.7) and ( 2.12 ) we obtain (see Appendix A) 

uj + uj 



2uj 



o J \v (w)r [tt 2 + z 2 (uj)} 



1 



uj — SI 

uj — UJ 

— : a u 

UJ + UJo 

1 2uj. 

UJ + SI UJ + UJ n 



+ z(uj)S(uj - Q) 



2uj n 



UJ + UJo 



v (SI) a u , 



—v* (SI) a u , 



(2.20) 

(2.21) 

(2.22) 
(2.23) 
(2.24) 



where 



4 



z(uj) 



lu 2 — lu 2 — 2lu H (lu) 
2oj \v (u>)\ 2 



and 



H (lu) = F (w) — G(lo)=V f \- ^ dSl - 

J uj — \l 

We can express a as a function of and in the following way 

du</> u A\, + / dujip^Alj. 



\y(m- 

u + n 



dSl. 



(2.25) 



(2.26) 



(2.27) 



Now taking, again, the commutators 



a, At 



and 



a, Au 



we obtain (f> u — a* and Lp u = — Xu- Substituting the 



expression (2.20) for A u in (2.27) the time evolution of the operator a can be easily written as 

dui 



where 



with 



a(t) 



\L (w)\ {A (uj) cos (cut) a-i[B(uj)d + C (lu) a f ] sin (tot)} 

r dfi , - r dsi , ~+ 

+ / — B 1 (Q;t)b n + J —B 2 (Sl-t)b\ l , 



A(uj) = 2oj, B(u) 



LU 2 + LU 2 



C(lu) 



LU Q LU a 



B-l (Si; t) = v (SI) {(lu + SI) [X (Si; t) + Z (SI) e" im ] - iY (+) (ft; t)} 
B 2 (0; t) = u* (O) { (u - SI) [X (0; i) + Z (Q) e j0t ] - iF ( _) (fi; *) } , 



X(Sl;t) 



du> 



Y { ±) (SI; t) = V J duo 



2|^MI 2 

LO 2 -SI 2 

io 2 -Sl 2 



LU COS (out) 



(lu 2 ± uu Sl) sin (u>i) , 



Z(Sl) 



\L(Sl)\< 
\v(Sl)f 



2 rr>2 



ft 2 -wi 



2w 



-H(Sl) 



(2.28) 

(2.29) 

(2.30) 
(2.31) 

(2.32) 
(2.33) 
(2.34) 



27TW |u (w)|^ 



[w 2 -w2-2w i/(w)] 2 + 2ttw |u(w)| 2 



(2.35) 



III. THE MODEL OF COORDINATE-COORDINATE COUPLING 

The expressions obtained for the evolution of the operator a(t), within or without the RWA, remained written in 
terms of the coupling function |u (w)| 2 . Therefore, the choice of the function \v (lu)\ 2 will determine the dynamics of 
the damped oscillator. We will choose the function \v (lu)\ by comparing the dissipation model corresponding to the 
Hamiltonian (2.4) to the one presented in pi that corresponds to the following Hamiltonian 



H 



P 

2M 



V(q) 



E 

i 



3_ 

2m 3 



(3.1) 



-^2,C } qjq + V R {qt), 
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where the counter-term Vr (q), which cancels the additional contribution to V (q) due to the coupling of the system 
to the reservoir, is given by 



C 



2 

3 ~ 2 f- 



2nijU)j 



(3.2) 



The spectral function J (lu) is defined by 



C 

j<w) = £ —i-s (w - wj) 



2 z — ' rrijUij 



* 9 (.<»)<% 

2 m,,,w 



where we have taken the limit of a continuous spectrum and used g (oj) from (EJ) . For ohmic dissipation 



J(w) = 



2M7W i/ w < O c 
if lu > Q c , 



(3.3) 



(3.4) 



where f2 c is a cutoff frequency much larger than the natural frequencies of the motion of the system of interest. But 
in our calculations we will conveniently use the Drude form 



J M 



2M7W 



(3.5) 



We are treating a damped harmonic oscillator so V (q) — l/2AIuj^q 2 . Applying the usual definitions of the operators 
a and bj, 



Mlo 
21% 



MiO n 



2h 



1i 



-Pi 



we can rewrite (3.1), initially without the inclusion of the counter-term Vr (q), as 

H — Tiuj a)a + huijbjbj 



(3.6) 



(3.7) 



2 V Muj 

(measuring the energy of the system from the energy of the vacuum). Now we can use the transformation ( |2.l|) in 
order to consider a continuous spectrum for the excitations of the reservoir. The second term on the RHS of (3.7) 
becomes 



Hrcs = h I Lubl.buidu) 



and its last term can be written in the following way: 

Now comparing (|3~7H3~9|) with Q we see that both Hamiltonians will be equivalent if we employ 



(3.8) 



(3.9) 



1 



2 y Mtu o m u L0 



Taking the square of (3.10) and comparing it with (3.3) we obtain 

1 J (to) 



v 1 M 



2tt Mu! ' 



(3.10) 



(3.H) 



G 



Adopting the Drude form (3.5) \v (w)| 2 is given by 

,2 7^ 



\v(u)\' 



1 



TTUJo (l + W 2 /fi 2 ) 



(3.12) 



which is defined only for uj > 0. 



Now that we have established the form of \v (uj) | corresponding to the Caldeira-Leggett model |2J , we can determine 



H (u) through (2.26). A simple calculation shows that H (uj) will be given by 

7^ c 1 



H(u) = - 



UJ (l+U, 2 /fi 2 ) 



(3.13) 



We can also diagonalize the Hamiltoni an ( 3.1) considering the inclusion of the counter-term Vr (q) (see Appendix 
A). The result is that all the equations ( 2. 20^2.35 ) will remain valid with the following substitution: whenever the 
function H (u>) appears it should be replaced by 



H R (uj) = H H 



2uj n 



where the frequency shift Auj 2 is defined as [E 



Auj 2 



1 N 

— Y 



c 2 



7=1 3 3 



\v(u))\ 7^ c 



L0 o 



(3.14) 



(3.15) 



Whenever a function appears with the sub-index r it means that we are considering the introduction of the counter- 
term. 

The spectral function ( [3.5| ) is appropriate to the description of the reservoir since we consider f2 c 3> uj ,j. So, in 
order to simplify and also obtain the exact function associated to the ohmic dissipation we will take the limit S^i c — > oo 
in the expression for \L (uj)\ 2 r . To do so, first we consider the renormalized function \L (uj)\ 2 r given by 

\L(.)\l= 2 ^ >M ' 2 r ^- (3.16) 



[uj 2 — uj 2 — 2uj Hr (ui)] 2 + 2-klJo \v (uj)\ 2 



Once 



we obtain 



lim Hr(u>) = 0, and lim |t>(a>)| z 

f2 c — >oo f2 c — >oo 



7<J 
TTUJ„ 



lim \L{uj)\ r 



(uj 2 -uj 2 ) 2 + (2juj) 2 ' 



Thus, we see that 



lim \L(uj)\ 2 r = M X " (uj), 



(3.17) 
(3.18) 
(3.19) 



where x (uj) is the imaginary part of the response function of a damped harmonic oscillator. 
In the limit 7 < w we can write 

2cj 7 



\L(uj)\] 



[2uj (uj - uj )] + (2w 7) 2 



1 



7 



2w Q (uj — ui ) + 7 2 



(3.20) 



that corresponds to a Lorentzian distribution of width 7. 

For the function \L (uj)\ 2 , without the renormalization, we have H (uj <C O c ) ~ —j£I c /uj and therefore for f2 c 3> u> , 7 
we obtain 

27a; 



\ L (u)\ = 



{uj 2 -uj 2 d + 2 1 n c ) 2 + (2 1 u) 2 ' 



In this case we should have uj 2 > 27S1 C , because, without the renormalization, we must have jl9[ 



uj 2 > \Auj 2 



(3.21) 
(3.22) 



for the diagonalization to be consistent. 
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IV. ANALYSIS OF THE EVOLUTION OF a (t) 



Now we can analyze in detail the ti me ev olution of the operator a associated with the system. We will analyze each 
term of the expression for a (t) in eq.( 2.28 ). We will be interested in the relation between the degree of dissipation in 
our system and the importance of each one of those terms. 

Initially we will analyze the coefficients associated to the operators a and a< . The fastest and most efficient way to 
understand the behavior of each one of them is through graphs. 

The graphs in Fig. 1(a) present the behavior of \L (uj)\ r for three values of 7: 71 = 0.1w o , 72 = uj and 73 = 10u> o . 
We see that for 71 = 0.lu) o , \L (w)\ R presents a narrow peak centered approximately about ui (we showed that in the 
limit 7 <C uj the function \L (uj)\ r tends to a Lorentzian centered at uj and with width 7). As 7 increases (72 = ui ) 
the function \L (uj)\ r broadens and becomes centered at progressively lower frequencies. For 7 still larger (73 = 10u> o ) 
\L (oj)\ r narrows again, but its peak is about very low frequencies. 

The graphs in Fig. 1(b) present the behavior of the functions A(uS), B(uS) and C (u>) that appear multiplying 
\L (<*))] R in the different terms of the expression for a (t). Simultaneously observing (l.a) and (l.b) we conclude that 
when 7«w„ the function C (ui) \L (uo)\ R has a negligible amplitude if compared to the functions A (uj) \L (uj)\ r and 
B (uj) \L (uj)\ r , because in this case \L (u>)\ R is very sharp and centered at uj whereas C (uj ) — 0. As 7 /u) increases, 
\L (oj)\ r has its peak broadened and moved away from uj . The function C (uj)\L (u)\ R becomes comparable to the 
others and in the limit 7 ^> ui , it is of the same order of B (uj) \L (uj)\ r whereas A (uj) \L (uj)\ r becomes very small. 



It remains to analyze the coefficients Bi tR and B 2 ,r (fyt) of b^ and b^, respectively, in the expression ( 2.28| ) 



for a(t). We know that in the limit 7 « 0) o the function |L(o;)| fl tends to a Loren tzian centered at lu and with 
width 7. Therefore, the function (uj — Q) Z R (f2) that appears in the expression (2.31) for B 2 ,r (O; t) is, in this limit, 
negligible if compared to the function (ui + Q) Z R ; (f2) in the expression ( 2.3C ) for Bi_r (fl;t). The evaluation of 
X R (fi; t) results in 



(n 2 



dt 



o: 2 



J 2 + 7 2 



sin (oj't) + 27 cos (uj't) 



-it 



27SI sin (nt) 

(n2-o;2) + (2 7 nr 



(4.1) 



for 7 < u) , where us 1 = y 'ui 2 — 7 2 . We see that in the limit 7 <C w Q the functio n X R (f2; t) will also be very sharply 



peaked around u) a . Therefore, the function (u — Cl) X R (Q; t) in the e xpression (2.31) for B 2 ,r (fi; t) is also negligible 
if compared to the function (oj + Q) X R (f2; t) in the expression (p0|) for B hR (tl; t). Similarly it can be shown that, 



in this limit, the function Yr_\ R (f2;t) is negligible in relation to the function Yr + \ R (p,;t). We conclude that in the 
limit 7 <C u>o the coefficient B 2 ,r (fi; t) is negligible in comparison to the coefficient Bi.r (f2;t). As the ratio r y/uj 
increases and the function \L (oj)\ r changes its shape, the coefficient B 2 r (£l;t) becomes comparable to B\ r (f2;t). 

So far we have analyzed the relevance of the terms associated to and 6q in the expression ( |2.28| ) for a (t) 
considering the inclusion of the counter-term Vr (q) in our model. We showed that these terms are negligible in the 
limit 7 <C u> , but become important as the dissipation increases and the function \L (oj)\ r becomes broader and is no 
more centered at ui Q . Now if we had not considered the inclusion of the counter-term in the interaction Hamiltonian, 
we would have |L(w)| 2 given by ( 3.21 ) instead of In this case, we see that the condition for the function 



\L (u>)\ to be centered very close to lj is that 27f2 c <C lo or 



U) fir 



(«!)• 



(4.2) 



Therefore, the condition j/ujo <C 1 would not be enough for us to ignore the terms associated to and b^ in the 
expression for a(t). These terms can only be neglected if the condition (4.2), which limits our system to a much 
weaker dissipation, is satisfied. 

We notice that a system subject to a weak dissipation (7 -C u> , in our case) does not guarantee that its frequency 
shift (Auj 2 = 27f2 c ) is also small. We will see later, in more detail, that for a system subject to very weak dissipation 
t he d amping coefficient 7 will be given by ir \v (w G )| and the frequency shift by H(ui )- Observing the expression 
(2.26) for H (uj) we clearly see that the relation between these functions depends on the form adopted for the function 
\v (uj)\ . Therefore, n \v (uj )\ <C uj does not guarantee that we will have H (uj ) <C u> (as we have seen to be the 
case for \v (u})\ 2 given by ( 3.12| )), although this can happen for some functions \v (uj)\ 2 . 
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V. REDUCTION TO THE MODEL WITH THE ROTATING WAVE APPROXIMATION 



Now let us consider the situation in which the following conditions are satisfied 



7r|w(w)| 2 «w , for u ~ ui Q 

H (uj) <C UJ , for UJ ~ LU 



(5.1) 
(5.2) 



Under these conditions the function \L (uj)\ 2 will be a function well peaked around ui a . Therefore we can ignore the 
terms associated with and 6^ in the expression for a (t) . Even the expressions for the coefficients of a and 6q can 
be approximated considering that \L (oj)\ 2 will only be appreciable, in this case, for uj ~ ui Q . We can write 



A ( w ) \L (uj)\ 2 ~ B (a;) |L H| 2 ~ 2w Q |i H| 2 , 



(5.3) 



and finally 



B. 



(i) 



vffllj duj2\L(uj)\ 2 V-^i 



a(t) 



\L(Sl)\ 2 ro 

-2cj 5- S2 — o> 

Mtt)| 2 



H(Sl)]< 



-int 



+ J dnv{n)S^J duj\a^v^- 



+- 



\v(Sl)\ 



[n - ui - h (si)] , 



(5.4) 



(5.5) 



where the function \a u \ comes from the approximation of |L(w)| considering (5.1 - 5.2) 



2 -^\L(uj)\ 2 
ir 



[uj — uj — H (uj)] 2 + ir\v(ui)\ 2 



(5.6) 



Now let us compare (5.5-5~^) with the expressions ( ^.14 ) and ( 2.19p , previously obtained in the RWA. The only 
difference between these expressions is given by the presence of H (uj) instead of F (uj). Once H (uj) — F (uj) — —G (uj) 
we would have, for uj ~ ui , H (uj) ~ F (uj) if G (uj) -C F (uj). There can be functions |u(w)| 2 that satisfy this 
r equir ement. However, most of the physically reasonable functions \v (uj)\ 2 do not; for example, if \v (uj)\ 2 is given by 
( 3.12 ) we have G (uj) /F (uj) ~ —1 for uj ~ uj . In this case, H (uj ) ~ 2F (uj ) yielding twice the frequency shift given 
by the model within the RWA |^2| . The same relation is found whenever \v (uj)\ 2 extends to frequencies much larger 
than uj with nonnegligible values, for then 



H (uj ) ~ -VP J dSl ~ 2F (lu ) . 



(5.7) 



This larger fre qu ency shift can be easily understood th roug h a perturbative analysis. Let us consider a system 
described by (1.1), (1.2) and having Hint within the RWA (1.4). It can be shown that, in second order, the perturbed 
levels of the oscillator remain equidistant with an apparent frequency uj + A rwa uj, where 16 



A 



RWA, 



(5.8) 



Taking the continuous limit and using (2.5) we see that this expression is nothing but F (w a ) which really represents 
the frequency shift in the weak dissipation limit. Now it is easy to show that if we consider Hi nt given by ( |l.3| ) 
without the RWA, we have in second order in the perturbation, 
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L0 o + LOj 



(5.9) 



This expression, in the continuum limit, is nothing but H(uj ). Therefore, we see that the substitution of F (u> ) by 
H (lu ) could already be foreseen by a simple perturbative theory. The same perturbative analysis can be used to 
understand why the counter-rotating term is not important in the calculation of the decay rate of the system in the 
weak dissipation limit. In first order, the decay rate of the system is given by the Fermi's golden rule for which only 
the terms of Hint that directly conserve energy in the transition are relevant. This is not done by the counter-rotating 
terms. In fact, it is only done by the rotating terms that create or destroy energy quanta such that Uj = lu . This is 
the reason for the dependence only on \v (to )\ that appears in the very weak dissipation calculations. 

In a model that takes the counter-term into account we automatically have Hr (uj) — and the expression (5.5) 
can be substituted by 



o(t) 



dVLv {n) 

\v(n)\ 2 



(Q - uj ) e 



(5.10) 



where 



2dJn , _ , ,,2 



— ^ — — jr 72 = \ a "\B ■ 

(u — u ) + n\v(w) ] 



(5.11) 



Therefore, the RWA leads us to the correct results, with regard to the decay rate of the system ^related to , 

if and only if the condition of weak dissipation ( |5.j[ ) is satisfied. Regarding the frequency shift (associated to F (w )), 
we see that its agreement with that given in the limit of weak dissipation, in a model without the counter-term, 
strongly depends on the function \v (w)| 2 adopted. For functions \v (uj)\ 2 that extend to frequencie s m uch larger than 



lu we have twice the shift foreseen in the RWA. Besides, it is also necessary that the condition (5.2) be satisfied in 
order to guarantee that this shift is much smaller than w (and we can neglect the terms in and Dq in the expression 
fora(t)). 

In the case of ohmic dissipation the conditions (5.1-5^) are reduced to 



7 < w 0) 



(5.12) 



once in this case 



7r|u(cj )| =7 and Hr (u>) = 0, 



in the limit fl c 



CO. 



VI. EVOLUTION OF A COHERENT STATE 



We showed that if our system satisfies the conditions of weak dissipati on (5.1) and small frequency shift (5.2) the 
evolution of the operator a (t) can be reduced to the expression given by ( |5.5| ). Now we will suppose that initially our 
system is in a coherent state |a) and that the reservoir is in the vacuum state |0) corresponding to a reservoir at zero 
temperature. In this case we have 



a(i)|a,0) = J 



duj 



\a,0) 



(6.1) 



Therefore, in this particular case, a coherent state stays as such during its evolution with eigenvalue a (t) given by 

a(t) = a { duj\& u \ 2 e~ luJt . (6.2) 
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We can also calculate the evolution of the operator Sq (t) of the reservoir. Then in the case of weak dissipation and 
small frequency shift we can show that the modes of the reservoir also evolve from the vacuum state to coherent states 
with eigenvalues given by 



(*) = « 



V / doj 



ft 



(6.3) 



ft- uj - h (ft) 2 _ int 



i«(«)r 



| an | e" 



(ft) . 



Still under the conditions (5.1-5.2) we can further approximate |a u | by 



7T \V {UJ )Y 



[u> - ll> - H {lo )] 2 + tt\v(uj )\ 2 



(6.4) 



and also extend the lower limit of the frequency integral in (6.2) to — oo introducing a negligible error. Then we have 

a (t) = ae -i["o+A u ]t e -*\v( Uo )\ 2 t^ where AuJ = H ( Wo j _ ( 6 5 ) 



In the case of ohmic dissipation with the inclusion of the counter-term we have 



a (t) = ae 



-iuj t -jt 



(6.6) 



Now it is clear that when (5.1-5.2) are not satisfied making the terms associated to the operators a) and b n in the 



expression ( 2.26 ) for a(t) no longer neglegible, |a,0) will not be an eigenstate of a(t) because \a) and |0) are not 
eigenstates of and &J>, respectively. Therefore, we see that an initial coherent state |a), interacting with a reservoir 
even at temperature T = 0, will not remain a coherent state during its decay unless we have a system subject to very 
weak dissipation. 

The previous works that emphasized the existence of dissipative coherent states Q-Q] , in models described by the 
Hint (1-3), were based on master equations obtained through a method that is appropriate only in the limi t of weak 
dissipation. However, we saw that in this limit the corresponding model ( jl.3| ) is reduced to the RWA model ( |l.4D that 
really preserves the coherent states. We believe that the implicit assumption of weak dissipation is the reason why 
these authors have obtained the dissipative coherent states. Our result agrees with the one presented in |l3| where it 
was shown that the model (|3.l|) presents the coherent states as the initial states of the system that produce the least 
amount of entropy as time evolves. 



VII. EVOLUTION OF THE CENTER OF A WAVE PACKET 



We can also study the evolution of t he o perator q associat ed to the position of the particle. Once the operators q 
and p are related to the operator a by (|3.6|), we obtain from (2.28) the following expression for q (f): 



where 



q{t) = Gs (q,p;t) + T n (qn,pa;t) , 



g s (q,p;t) =q^-£(t) + j^£(t). 



— V (ft) 



too ft 
Muj 



dt 



W R (ft, t) + Z R (ft) cos (fit) 



?n 



[VLW R (ft, t) + Z R (ft) sin (fti)]} 



/d 
^|LM| 2 fl sinM), 

w R (n,t) =■ 



Pn 
mnft ' 



da; ' v /,'f sm (tot), 



(7.1) 
(7.2) 

(7.3) 
(7.4) 

(7.5) 
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with Z (Q) defined in (\U§. 

Now we suppose that the initial density operator of our global system can be written in the factorizable form 

Pt = Ps® pn, (7.6) 

where p$ and pn are, respectively, the density operators of the system and reservoir when they are isolated. Then we 
have 



(q (t)) = Tr s [Q s (q,p; t) p s ] + Tr n [T n (qn,PiY, t) pn] 
= Qs ({q) s > (P)s > t) ~\~ -Ftz {{qn) n , {Pn) n ] t) . 
Assuming that the initial state of the reservoir is such that 

%)n = (Pi)n = °> 

which in the continuum limit corresponds to (<?n) K = (Po)tc = 0, we obtain the following expression for (q (t)}: 



(q{t)) = (q) s j t jC(t) + ^fC(t), 



where 



^sin (w't)e-T*, 
C(t) = { t.e-^, 

-tit _|_ 



72—71 



71—72 



for 
for 

" 72t , for 



7 < w , 
7 = w , 
7 > Wo, 



(7.7) 
(7.8) 

(7.9) 
(7.10) 



with oj' — y/u>g — 7 2 and 7^2 = 7 ± a/7 2 — The expression (7.9) was also obtained by Grabert and collaborators 
P), by the method of functional integration. They affirmed that it would correspond to the classical trajectory of a 
damped harmonic oscillator. However, it is easy to see that this is not true. If the initial state of th e sy stem presents 
an initial average momentum (p) s — p and an initial average position (q) s = q a , then according to ( [7.9] ) (q (<)) would 
evolve as 



/ 7 

(q (<)) = q cos (uj't) sin {ui't) 



- 7 t 



(7.11) 



Po 

Mlo 



- sin (u> t) e 



-it 



for 7 < lo . However the classical trajectory is known to be 



7 

1 Wdas = Qo cos (u't) + — sin (ut) 



-ft 



(7.12) 



+ ^ sin ( w,t ) e 



-■yt 



Muj 



Thus, we see that there is a phase difference between ( 7.11 ) and ( 7.12 ) if the oscillator has an initial displacement q . 
Let us now suppose that the initial state of the reservoir is such that 



(?) 



(P. 



0. 



(7.13) 



We can write the expression (7.3) in the discrete limit, replace (7^8) by ( 7.13 ) and return to the continuum limit. 
Then we obtain (see Appendix B) the following expression for (q (t)): 



(«(*)> = <& 



-£(t) + 2 1 £(t) 



(P)s 
M 



C(t). 



(7.14) 



In this case if the in itial state of the system presents an initial average momentum (p) s = p and an initial average 
position (q) s = q a , (7.14) becomes 



7 

(q (<)) = q a cos (uj't) H sin {ui't) 



,—yt 



(7.15) 



for 7 < ui , that corresponds to the correct classical trajectory. 

Thus, we see that the classical evolution is not obtained with the initial condition ( |7.8| ) but with the initial condition 
(ff7l3| ). We can understand why this happens through the classical analysis of the model ( |3.5[ ) presented in the next 
section. 
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VIII. CLASSICAL ANALYSIS AND DISCUSSION 



In this section we will accomplish a classical analysis of the model used . Our obje ctive is to obtain a physical 
intuition on the effec t tha t causes the difference between the equations ( 7.11 ) and ( 7.15 ) and then on the meaning of 
the initial condition ( [7.13j ). This procedure can be justified by the equivalence of the classical and quantum dynamics 
of this model @. 

The Hamiltonian ([O]) can be written as |2lj| : 



2m, 2 



Qj 



C 3 

— —Q 

m J cu J 



(8.1) 



The equations of motion of this system are given by 

Mq(t) + V'(q) =J2 C J 

3 

rrijiij (t) + rrijU -qj (t) — C 3 q (t) 



Qj (*) 



C, 



!«(*) 



If qj (0) and qj (0) are the initial conditions the solution of the homogeneous part of fl8.3| ) will be 

q f (t) = q } (0) cos fat) + sin fat) . 

Ulj 



(8.2) 
(8.3) 

(8.4) 



Th e pa rticular solution, considering the presence of the force Cjq (t), can be obtained by taking the Fourier transform 
of (8.3). Then we have 



tI{t) = ^k$ldt'q{t')sHu j {t-t>)] 



C 



ffjT {q (t) - q (0) cos (ujt) 

-f*dt>q(t>)co S [cUj(t-t')]} 



•5) 



Using the definition of the spectral function J (uj) (|3 . 3|) and (3.4) it can be shown that in the limit O c — ► oo we have 



E 



c 2 



m uj j j 



dt' cos [uj (t ~ t')] q (*') = 2Mjq (t) . 



(8.6) 



Therefore, the general solution of (S.3), qj (t) = q^ (i) + q? (t), when substituted in (8.2) results in the following 
Langevin equation 



where 



Mq(t) + V (q) + 2M iq (t) = F(t). 



F w = E c & (°) cos to*) + J2 —to (°) sin to*) 



•7) 



is the fluctuating force and we have redefined the position of the oscillators of the bath |2c 



Qj (0) = Qj (0) 



9(0). 



(8.8) 



Supposing that the bath is initially in thermodynamic equilibrium in relation to the coordinates qj (0) we have, in 
the classical limit, 



(&(0)) =(q j (0)) = (q j (0)q f (0))=0. 
% (0)^(0)) = ^%>, (qj (0)^(0)) " 



m 



(8.9) 
(8.10) 
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The physical meaning of this initial con dition wr itten in terms of the relative coordinates qj has ahead been analized 
by Zwanzig |^4| some time ago. Using (B.S- B.10 ) and after some algebraic manipulations it is shown that (F (t)) = 
and (F (t) F (£')) ~ AM'jkTS (t — t') which correspond to the expressions that characterize the Brownian motion. 
On the other hand if we had adopted the initial condition 



(0)) = (<i 3 (0)) - ^ (ft- (0)) = 0, 



(8.11) 



we would have 



(F(t)) = -3(0)^^ cob (w 3 -t) 

^ ' TY1. -It). 



-4M73 (0) 



cos = -4M~/q (0) 5 (i) . 



.12) 



where we have used (3.3), (3.4) and taken the limit fl c — > 00. Therefore, we would not have (F (i)) = 0, but the 
presence of a delta force at t = 0. Physically what happens is th at if the oscillators of the bath are not "appropriately" 
distributed around the particle (as in the initial condition ( |S.ll| )), when it is inserted in the bath, these oscillators will 
"pull" the particle until they reach this "appropriate" distribution. This force will act on the particle during a time 
interval of the order 1 /Q c . Therefore, in the limit Q c — > 00 we will have a delta force that will ca use a phase difference 
in the evolution of the system. This phase difference is the difference betwe en (|7.1l| ) and ( 7.12|) which is corrected in 
( |T. 1 5|) b y the adoption of the initial condition (7.13) (quantum anal ogue of (|3.9|)) inst ead of ( |7.8| ) (quantum analogue 
of QB.11|)). As far as we know, the need to use the initial condition ( |7.13| ) in place of in the quantum treatment 
of this model has not been noticed in previous works. In Ref. 0] the authors make some approximations which are 
equivalent to regarding the initial time as t = + (t ~ l/f2 c ). So the initial conditions are established at this instant 
although the coupling between particle and bath is switched on at t = and gives rise to a delta type force at this 
instant. The inclusion of t = in propagator methods must be accompanied by the above-mentioned modification 
of the factorizable initial condition. However, it must be emphasized that we are not addressing here the question of 
the generalized initial condition p,^0|. Actually the point we have raised is clearly responsible for the disagreement 
between (q (t)) found in these references. In Ref. H the authors re prod uced the dephased (q (i)) (c.f. eq. (7.S) above) 
whereas in Ref. |L(| this time evolution is the correct one as in ( ffj5l ). The origin of the discrepancy is the use of 
t = or t = + as the initial instant together with the factorizable initial condition. 

We would like to take advantage of this opportunity to correct a mistake that was made in Ref. jl(| of which one 
of us is co-author. The referred article considers an initial condition of the system when the bath of oscillators meets 
thermodynamical equilibrium with the particle at the position it is placed in the bath. In this case one obtains mean 
values of the position (q(t)) and momentum (p(t)) which depend on the temperature of the reservoir and that do 
not exactly coincide with their classical counterparts. This disagreement was justified within a classical analysis of 
the model. In this analysis it was affirmed that the classical initial condition equivalent to the proposed quantum 
initial state, that corresponds exactly to (qj (0)) = 0, would imply in a classical solution of the model different from 
the trajectories of a damped harmonic oscillator. We saw in the present work that this is not true and therefore this 
argument can not be used. We believe that the origin of the disagreement when adopting a non-factorizable initial 
condition is the impossibility to desc ribe the evolution of the system through an independent sum of functions of the 
system and reservoir variables as in (7.7). The quantum effects of the correlation between the variables of the system 
and reservoir prevent a direct comparison of the quantum mean values with the values obtained through the classical 
analysis of the model. Accordingly, it can be shown that the discrepancy vanishes in the classical limit [kT 3> fuv ). 

After we had made the above analysis we became aware that in previous works pTLOl the authors had also noticed 
the existence of initials kicks and jolt s in this system when the initial condition ( |7.6| , [7.8[ ) is used. In both works the 
existence of an initial kick, given by ( 8.12 ) in the limit limite il c — > 00, is noticed (eqs. (3.2) and (45), respectively). 
However, the existence of this initial transient is considered as a character istic o f the model to be taken into account. 
In our analysis we see that although the existenc e of the kick given by ( 8.1 2| ) is a real characteristic of the model 
when it is subject to the initial condition (7.6,7.8) it is an undesirable feature that should be corre cted. Fortunately, 
this correction can be made even with an improved factorizable initial condition, that is, (7.6, 7.13] ). 

The authors of |[l| and [Q have recognized that the presence of initial jolts, in their master equation coefficients, 
generates certain non-physical effects and so they suppose that they are due the adoption of a factorizable initial 
condition. In jl4| t he evolution of the system is analyzed for a non-factorizable initial condition, similar to the one 
used in || and |10|], in which the initial position of the particle is defined by a measurement process in a state in 
thermodynamic equilibrium with the bath. However, the initial jolts in the time scale l/f2 c still persists. We believe 



that this happens because this initial condition does not satisfy (7.13) when the initial mean values of the position of 
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the particle and the oscillators in the bath are calculated. Thus, in this aspect, it is less general than the improved 
factorizable initial condition that we considered. 

Actually the initial jolt, at least the one they attribute to the decoherence process (however see discussion below), 
does not appear in the more general initial condition later adopted in [|15f . It is an initial condition prepared by 



a dynamic process in a finite time t p . In this case we can consider that the condition (7.13) will be satisfied since 
t p 3> l/f2 c - Indeed in this situation (t p 3> l/f 2 c ), it was shown that the initial jolt does not appear. 

Thus, we believe that the initial condition ( |7.13 ) is enough to eliminate most of or maybe all the initial transients 



that would appear in this system in the characteristic time scale l/f2 c . However, another initial transient in this system 
is also known. In |2j| it was shown that for a factorizable initial state in the high temperature limit (kT 3> 7if2 c ) 
of the master equation it presents an initial transient within the time scale of the internal decoherence of the initial 
wave packet. If applied to times shorter than this it can lead to nonsensical results. We believe that this pathology 
can only be really corrected with the adoption of non-factorizable initial conditions. 

IX. CONCLUSION 

In this paper we have applied the Fano diagonalization procedure to two Hamiltonians commonly used as models 
for dissipative systems in quantum optics and in condensed matter systems; the rotating wave and the coordinate- 
coordinate coupling models, respectively. 

By exactly diagonalizing these two models we have succeeded in showing how the RWA turns out to be the extremely 
underdamped limit of the more general coordinate-coordinate coupling model. We have also been able to analyze the 
role played by the counter-term in this limiting procedure from the latter to the RWA. We have sho wn through the 



evaluation of the destruction operator a (t) of the system that the RWA is a good approximation for (1.1 ) if and only 



if the conditions (|5.l|) and (5.2) are satisfied. For certain choices of |w(w)| 2 , we have H (u> ) « 7r|w(w )| z and the 



2 



fulfillment of ( |5.l[ ) automatically implies (5.2). However, for other choices, we can have H (u> ) 3> tt \ v (w q )| 2 and ( |5.2| ) 
limits the validity of the approximation. Once these conditions are satisfied, we have shown that the time evolution 
of the system is identical to that determined within the RWA, with the exception of the frequency shift. We have 
found that this shift will be given by H (u ) instead of F (u ). As we have shown these functions usually have the 
same order of magnitude, but they are not identical. For functions \v (uj)\ that extend to frequencies much larger 
than oj we have H (ui ) ~ 2F (w )- 

The comparison of the Hamiltonian ( [2.4| ) with the Hamiltonian of the coordinate-coordinate coupling model es- 
tablished the relation ( |3.1l| ) between the spectral function J (ui) of this model and the coupling function \v (w)| 2 . In 
the case of ohmic dissipation and considering the inclusion of the counter-term, we find that Hr (lu) = in the limit 
17 c — ► oo. Then the only condition required for the RWA to be valid is 

7<^o- (9-1) 

As an application of this method, we have studied the existence of dissipative coherent states and concluded that 
it can only exist within the RWA and when thermal fluctuations are neglegible. When these conditions are not met, 
the initial state will in the long run become a statistical mixture. 

Finally, we have also addressed the question of the discrepancies in the time evolution of the observables of the 
system that arise when the factorizable initial conditions are not properly accounte d for. We have sho wn h ow to deal 



with this problem by using the appropriate improved factorizable initial condition (7.13) rather than (7.8). 
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APPENDIX A: DIAGONALIZATION WITHOUT THE RWA 



Here, the procedure used in the diagonalization of the Hamiltonian (2.4) will be presented. We want to find the 
operator that allows us to writ e ( | 2.4| ) in the diagonal form. We write A u in its general form ( 2.20D and then we 
impose the commutation relation (|2.7jj 
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A.„H 



(Al) 



Replacing (2.4) and ( 2.20 ) in (Al) and taking the commutators of the expression obtained with a 1 , a, bn and b^, 
we have, respectively, 

(A2) 

(A3) 

(A4) 
(A5) 



uja^ = ujoa^ + J \j3a, t nv* (Cl) - a^^v (Cl)] dCl, 
luxuj = -UoXuj + / \p u ,aV* (Cl) - &u,nv (ft)} dCl, 



= (at u - Xuj) v (Cl) + Cl/3 Uj n, 
u<JuM = (a u - Xu>) v* (Cl) - Cla^fl- 



Subtracting (A2) from (A3), we have 



UJ + UJo 



Replacing (A6) in (A4) we obtain 



V 



1 



uj — Cl 



+ z(uj)5(uj - Cl) 



2w 



UJ + UJo 



—v (Cl) a. 



where z (ui) is a function to be determined. Similarly, substituting (A6) in (A5), we have 



UJ + Cl UJ + UJ„ 



2^o * /rA 
-v (J2)a w - 



Now, substituting (A7-A.8) in (A.2) we obtain z (uj) given by (2.25J 



It remains to determine a u . For this we impose the condition ( |2.12 ) which results in 

a u a% + J rffi/3 w ,n/3S,n ~ X^xS 

- / dCla^nal^ = 6 (uj - Cj) . 



Using (|A6|) and (|AS|) we obtain, respectively, 



2uj (uj + uj) 

(UJ + UJo) (Cj + UJ ) 



and 



(2uj ) 2 G(uj)-G(uj) 

(uj + UJ ) (UJ + UJ ) UJ — UJ 



were G (uj) is given by ( 2.26 ). Now, using (A7), as well as the property 
V V 1 ( V V 



UJ — UJ' UJ — UJ' UJ — UJ \ UJ — UJ' UJ — UJ 



- ) + tt 2 5 (uj - Cj) 8 
J 



uj' --(uj + uj) 



we obtain 



dCl/3 w ,Q0~n = 



(2uj 



(uj + UJ ) (0 + UJ ) (uj — UJ 



Cj 2 -uj 2 
2uj 



+ G(uj)-G(uj) 



W 2 + z 2 ■ 



(uj)] \v(oj)\ 2 S(oj-oj)} 
Then substituting flAlOHAllD and (|A13| ) in jA9| ), we have 

[tt 2 + z 2 (uj)] S (uj - Cj) = 5 (uj - Cj) 



* (2w D ) \v(w)\' r _ 2 .2 



' (UJ + UJ ) (UJ + UJ ) 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 



(A12) 



(A13) 



(A14) 
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and, therefore, we should have \a u \ given by (2.21). 

In the calculations presented above we supposed that \v (io)\ is a continuous function and such that \v (0)| = 0. In 
this way we guarantee that J °° dSlf (SI) \v 5 (Si — uj) — f (uj) \v (uj)\ for any nonsingular function / (uj) within 
the whole interval (0,oo). 

We can also diagonalize the Hamiltonian (^1) considering the introduction of the counter-term Vr (<?)■ Rewriting 
it in terms of the operators a and bj, defined in (3J3), we have 



H = fko a'a + h 



4w 



— (a + a T ) 



V hwMj ~ % < U- (a + at) V (bj + b}) 



(A15) 



Writing ( A15 ) in th e co ntinuum limit and following the same procedure as adopted above, we will see that the 
equations flA2| ) and (A3) will be substituted now by the equations 



u>a u 



Am 2 

2uj 



Am 2 
2u) 



+ / [Pu,nv* (A) " <r u jav (SI)] dSl, 



2U! 



2uj 



(Xuj + [Puj.nv* (SI) - <j w ,nv (SI)] dSl, 



(A16) 
(A17) 



respectively. The equations (A4) and (A5) will stay the same. Thus, it can be easily shown that all the other previous 
equat ions will not change with the only difference that the function H (uj) should be substituted by Hr (uj) given in 



APPENDIX B: CALCULATION OF Ttz , (pj) K ;t) 

The expression ( |7.3| ) for T-ji (qn,pa;t) can be written as 

Pn 



J(Sl-t) q n + JC(Sl;t) 



mnSl 



(Bl) 



where the expressions for J (S I; t) and K, (Sl;t) are obtained by direct comparison between (Bl) and (7.3). Now we 
can substitute the expression ( 3.10 ) for v (SI) in (Bl) and write the expression obtained in the discrete limit 



Fn(to,Pn;t) =-i-V^i J (Slj;t) Jgjsi/) [ dSlq, 



mn.Slj 



g(Slj) 



dSlps 



Recalling the relation (2.1) between the discrete and continuous operators, we obtain 



?n (Qj,Pj;t) 



M ^ IT 

3 



J (Slj\ t) qj + K, (Slj ; t) 



Pj 



mnSlj 



Employing the initial condition (7.13), we have 

Fizdq^RAP^n^) 



Mir ^ 

nit) (q) s , 



Mn ^-r 1 m.jSl 2 
j j 



C 2 

; J(Slf,t)(q). 



with 



n (t) = -auj 



dSl \v(Sl)\' 

~k sT~ 



J(Sl;t). 



(B2) 



(B3) 



(B4) 



(B5) 



where we used again the relation (3.1C). Writing Ti (t) as 
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H (t) = h (t) + h (t) , (B6) 

we have 



h («) = -4o; / ^^L Zr (fi) cos (fif) . (B8) 

J TT Q 

The calculation of I\ (t) is a somewhat lengthy but straighforward calculation and results in I\ (t) = 0. So all that is 
left is 

r do, lu 2 - n 2 

H (t) = I 2 (t) = 4 7 / 2 cos (Sit) . (B9) 

J TT (n 2 - u 2 f + (2jn) 2 

The evaluation of this last integral can also be accomplished by the method of residues and yields 

H (t) = 2 1 C (t) , (BIO) 

for t > 0. 
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FIG. 1. (a) Graph of |L(w)|^ for different rations "f/oJo- (b) Graph of the functions A(lu), B (oj) and C (uj) that appear 
multiplying |L (u)|^ in the different terms of the expression for a (t). 
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